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Abstract 

(Et4N)3Fe2Fg exhibits a butterfly-shaped hysteresis below 5 K when the magnetic 
field is parallel to the threefold axis, in accordance with a very slow magnetiza- 
tion relaxation in the timescale of minutes. This is attributed to an energy barrier 
A = 2.40 K resulting from the S = 5 dimer ground state of [Fe2Fg] 3_ and a negative 
axial anisotropy. The relaxation partly occurs via thermally assisted quantum tun- 
neling. These features of a single-molecule magnet are observable at temperatures 
comparable to the barrier height, due to an extremely inefficient energy exchange 
between the spin system and the phonons. The butterfly shape of the hysteresis 
arises from a phonon avalanche effect. 
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1 Introduction 



The discovery of macroscopic quantum spin tunneling and molecular mag- 
netic hysteresis at cryogenic temperatures in the spin cluster compounds 
[Mn 12 12 (OAc) 16 (H 2 0) 4 ] -2HOAc-4H 4 (Mn 12 ),[l-4] [Fe 8 2 (OH) 12 (tacn) 6 ] 8 + 
(tacn = 1,4,7-triazacyclononane) (Fes), [5] and various Mn 4 clusters [Mn 4 3 X 
(OAc) 3 (dbm)3] (dmb - = dibenzoylmethanate, X = F, CI, Br, OAc) (Mn 4 )[6,7] 
has revitalized the field of molecular magnetism. Chemists have joined forces 
with physicists in an attempt to fully characterize and understand the new 
phenomena and to develop strategies to design and prepare other members of 
this new family of single molecule magnets (SMMs). So far a number of fur- 
ther SMMs have been discovered, mostly by means of ac susceptibility, among 
them several manganese clusters, [8-11] as well as tetranuclear vanadium[12] 
and iron [13] complexes. All these compounds contain clusters with at least 
four paramagnetic centers. 

Here we report low-temperature magnetic properties of (Et 4 N) 3 Fe 2 F 9 (Fe 2 ), 
which clearly show effects of slow relaxation and quantum tunneling. To our 
knowledge the [Fe 2 F 9 ] 3 ~ dimer is the smallest exchange-coupled unit so far to 
exhibit this type of phenomenon. 

2 Experimental 

(Et 4 N) 3 Fe 2 F 9 was prepared as small needle-shaped crystals according to ref. [14] 
The product was characterized by X-ray powder diffraction using the struc- 
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tural information from ref. [15]. It crystallizes in the space group P63/UL with 
the trigonal axes of the [Fe 2 F 9 ] 3 ~ dimers parallel to the hexagonal c axis, see 
the inset of Fig. 1. Due to their water sensitivity, single crystals of Fe 2 were 
sealed in glass capillaries for the magnetic measurements. The crystals were 
perfectly transparent, indicating that no hydrolysis occurred. 

Magnetic measurements were performed using a Superconducting Quantum 
Interference Device (SQUID) magnetometer (Quantum Design MPMS-XL-5 
with a 5 T (50 kG) magnet), with the magnetic field H parallel and perpen- 
dicular to the crystallographic c axis. 



3 Results 

The single-crystal magnetic susceptibility is shown for H \\ c and H _L c 
orientations in Fig. 1. Below 25 K it is highly anisotropic, and was fitted using 
the spin Hamiltonian 

H = -S 2 ) + DJ2 \sl - l(Si(Si + 1)1 + E MBS* ■ H (1) 



i=l,2 



i=l,2 



with J = -2.23 K, D = -0.215 K, and g = 2.00. The exchange splitting is 
ferromagnetic with an S = 5 dimer ground state, separated by 11.2 K from the 
next higher S = 4 state. The single-ion axial anisotropy parameter D can be 
correlated with the zero-field-splitting (ZFS) parameter of the dimer ground 
state by D s=5 = 4/9 ■ D = -0.096 K.[16] 

Magnetization measurements with H _L c showed the expected behavior for 
the above Hamiltonian and parameters, see the inset of Fig. 2. The results 
of magnetization measurements with H || c (H z ) were unusual. Fig. 2 shows 
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1.8 K results obtained with two different sweeping rates T. They exhibit a pe- 
culiar butterfly-shaped hysteresis with no remnant magnetization at H z = 0. 
Similar behavior has been observed in Mn 12 at 3.25 K, i.e. just above the 
blocking temperature. [17] The breadth of the wings increases with decreasing 
temperature and with increasing T. This is direct evidence for a slow relaxation 
process at 1.8 K. The effect is strongly T dependent, and we could observe it 
up to 5 K. For an infinitely slow T the butterfly hysteresis would disappear 
and the adiabatic curve (full line in Fig. 2) would be obtained. In order to 
quantitatively determine the relaxation dynamics, we performed magnetiza- 
tion relaxation experiments at fixed magnetic field H z and temperature T, 
after saturation at 30 kG and subsequent quick reduction (r ~ 150 G/s) to 
the indicated H z value. Typical decay curves are depicted in Fig. 3a for some 
selected values of H z and T. The relaxation depends on both. At 1.8 K and 
H z = 3.3 kG the relaxation is so slow that the equilibrium magnetization has 
not been reached after 40 minutes. For H z = at 1.8 K, on the other hand, 
the decay is about 4 times faster. At all fields the relaxation rate increases 
with temperature, as shown for T = 1.8 K and 3.0 K at H z = in the lower 
panel of Fig. 3a. Relaxation times r were obtained from single-exponential fits 
to the relaxation curves for various H z and T values. The field dependence of 
r at 1.8 K, shown in Fig. 3b, exhibits 3 distinctive dips at H z = 0, 1.1, and 
2.7 kG. The temperature dependence of r for a given value of H z was fitted 
with the Arrhenius law 

t = r ■ exp{AE/kT) . (2) 

At H z = the obtained parameter values are r = 70 s and AE = 2.2 K, 
corresponding to r = 238 s at 1.8 K. 



4 



Fig. 4 shows magnetization data obtained as follows: After saturation at +30 
kG the field was quickly (r 150 G/s) reduced to the given H z values, 
including negative ones. The magnetization was measured within 50 s after 
reaching H z . During such a fast sweep from +30 kG to —2.5 kG, e.g., the 
magnetization switches instantaneously from positive to negative saturation 
when crossing H z = 0. 

4 Analysis and Discussion 

Except for the butterfly hysteresis, the low-temperature magnetization curves 
in Fig. 2 are perfectly reproduced with the parameter values derived from 
the DC magnetic susceptibility data in the temperature range 1.8 K - 300 
K. Furthermore, the shortest Fe-Fe distances between neighboring dimers are 
at least 8 A, [15] which makes interdimer interactions extremely inefficient. 
Therefore the observed slow relaxation and butterfly hysteresis must be a 
genuine property of the [Fe 2 F 9 ] 3 ~ dimer. 

Our interpretation is based on a comparison with the observed behavior of 
single-molecule magnets such as Mni 2 [l,3], as well as the formalism developed 
for their interpretation. The appropriate picture representing the ZFS of the 

5 = 5 dimer ground state with the negative D s=5 is shown in Fig. 5. The 
barrier height in Fig. 5 at H z = is A = (\D S=5 \)S 2 = 2 AO K. This is 
slightly larger than the kinetic activation energy AE = 2.2 K determined from 
the magnetization relaxation. This difference is small, an indication that, if 
tunneling processes do contribute to the relaxation at the lowest temperatures, 
this contribution is not a dominant one. In the C^h dimer symmetry, the full 
anisotropy Hamiltonian includes the higher-order terms -B4O4, BqOq, and 
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S|6|.[18] The BlOl term mixes wavefunctions with AM S = ±6. This allows 
for resonant tunneling at H z = between the Ms levels — 3/ + 3 and at 
applied fields when the M s = —2/ + 4 and — 1/ + 5 levels cross, leading to 
significantly lower values of r at these fields and hence to dips in a r vs H z 
plot. Therefore the observed dips in Fig. 3b are direct evidence for tunneling 
processes in Fe 2 . The dips are located at H z = 0, 1.1, and 2.7 kG instead of 
0, 1.4, and 2.85 kG, respectively, calculated without the higher-order terms. 
This indicates nonzero values for the parameters B® and B®. In Fe 2 , the drop 
of r due to resonant tunneling compared to pure thermal relaxation at non- 
resonant fields is only about 25-75%, see Fig. 3b. In contrast, r decreases by 
1-2 orders of magnitude in Mni 2 when resonant tunneling is possible. [3] This 
is in line with our conclusion from the kinetic energy barrier, that although 
they occur, quantum tunneling processes are not dominant in Fe 2 . Instead, 
the purely thermal relaxation is dominant at all temperatures and field values. 
This distinctly different behavior from other SMMs is due to the combination 
of two facts: i) the thermal barrier in Fe 2 does not exceed the barrier for 
tunneling processes by more than a factor of 1.6; and ii) the tunneling is 
phonon assisted and thus slowed down by an extremely slow energy exchange 
between the spin system and the phonon bath. This latter effect will now be 
analyzed. 

The prefactor r = 70 s obtained from the Arrhenius analysis of the relaxation 
dynamics at H z = is unusually large, 7 — 11 orders of magnitude larger 
than in Mni 2 [l-4] and other SMMs [5,6,8-13]. It is much larger than the 
timescale of ac magnetic susceptibility measurements, which are usually used 
to characterize and quantify the relaxation behavior in SMMs. Only due to this 
extremely high tq value are we able to observe a slow relaxation phenomenon 
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below 5 K with a barrier height of only A = 2.40 K. In the SMMs reported 
so far, for comparison, A ranges from 7.6 K in [Mn403F(OAc)3(dbm)3][7] 
to 79 K in (PPh 4 )[Mn 12 12 (0 2 CEt) 16 (H 2 0) 4 ][9]. As indicated in Fig. 5, the 
magnetization in Fe 2 can relax both by thermal activation over the energy 
barrier (full arrows) and by tunneling between the Ms = — 3/+3 levels (dashed 
arrows). Both mechanisms rely on an interaction between the spin system 
and the crystal lattice degrees of freedom, i.e. on spin-phonon coupling. In 
this process, energy quanta corresponding to the energy differences between 
adjacent spin levels are exchanged between the spin system and the phonon 
bath. The extremely large tq in Fe2 results from an extreme inefficiency of 
these processes below 5 K, and this can be ascribed to the following two 
factors: a small spin-phonon coupling parameter g and a small density of 
phonon states in resonance with the energy spacings in Fig. 5. Quantitatively, 
the rate constant of a AM S = ±1 transition induced by phonon absorption 
or emission can be expressed as: [19] 

u/ _ 9qS±i (£m s ±i - £m s ) 3 /„x 

W Ms ±1,Ms ~ 487rpv , h 4 e (e Ms±1 -e Ms )/kT _ j > V) 

where S±\ = (S =F Ms)(S ± Ms+i)(2Ms ± l) 2 . The number of phonons in- 
volved in the process depends on their energy {sm s ±i — £ m s ) an d the energy 
distribution of the density of states, expressed by the sound velocity v in 
eq. 3, the mass density p = 1.36 g/cm 3 [15] and the Bose population factor 
( e (eM S ±i-e Ms )/kT _ ^^ 19 ] We note that ( £Ms±1 _ £Ms } i s about 1 order of 

magnitude smaller in Fe 2 than in Mni 2 and other SMMs. In all the SMMs 
reported so far, [1-13] the clusters are relatively large in size and contain or- 
ganic ligands such as acetate. This gives rise to a large number of vibrational 
and rotational low-energy internal modes, leading to a high density of states 
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in the energy range most relevant for spin-phonon coupling. This is in sharp 
contrast to Fe 2 , in which the [Fe 2 F 9 ] 3 ~ unit consists of atomic ligands with 
stiff Fe-F bonds, reflected in the absence of low-energy vibrational or rota- 
tional modes. Compared to the other SMMs, the stiffness of the [Fe 2 F 9 ] 3 ~ unit 
results in a strongly increased average sound velocity v in the crystal. Thus in 
Fe 2 the density of phonon states p p h oc l/t> 3 [20] is extremely low in the rele- 
vant energy range. The magnitude of the spin-phonon coupling parameter g 
is governed by the strengths of spin-orbit and electron-phonon (orbit-lattice) 
couplings. [21,22] The latter is expected to be smaller in Fe 2 compared to other 
SMMs due to the great stiffness of the [Fe 2 F 9 ] 3 ~ unit. As spin-orbit coupling - 
in terms of the ZFS parameter D s for the ground spin state - is also smaller in 
Fe 2 compared to other SMMs, [1-13] we can expect a smaller spin-phonon cou- 
pling constant go in Fe 2 . Thus the efficiency of spin-phonon coupling is greatly 
reduced in Fe 2 compared to other SMMs, and the smaller {sm s ±i — ^m s ) 3 j the 
larger v 5 , and the smaller g^ factors all contribute to the drastic reduction of 
the transition rate constant W in Fe 2 . As a consequence, the prefactor r in 
the Arrhenius expression is increased by 7-11 orders of magnitude. 

The magnetization data shown in Fig. 4 obtained with T m 150 G/s exhibit 
a very fast reversal of the magnetization around H z — 0. We ascribe this to a 
so-called phonon avalanche, [20] arising from a resonance of phonon absorption 
on one side of the barrier with phonon emission on the other side. A similar 
phenomenon has been observed in Mn 12 at H z ^ 0, induced by heat pulses[17] 
or extremely high sweeping rates. [23] The butterfly shape of the hysteresis in 
Fe 2 is very similar to that reported for K 6 [Vi 5 As 6 042(H 2 0)] • 8H 2 0,[24] but 
in the latter system it has been attributed to a phonon bottleneck arising from 
dissipative spin reversal in the absence of an energy barrier. Our results show 
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that a butterfly-shaped hysteresis can also occur in the presence of a barrier. 

Due to the phonon avalanche in Fe 2 , the reversal of the magnetization at 
H z = is a much faster process than its relaxation to equilibrium. The smaller 
the sweeping rate T in the magnetization experiment, the closer is the system 
to equilibrium, as is nicely demonstrated by the decrease of the butterfly hys- 
teresis with T in Fig. 2. After a fast (T m 150 G/s) reduction of H z to 0, the 
system is therefore still out of equilibrium at the beginning of the relaxation 
measurement. Since its reversal is so fast, the magnetization then overshoots, 
i.e. it becomes slightly negative, before it relaxes to its adiabatic value. There- 
fore there is an observable negative magnetization relaxation at H z = 0, see 
Fig. 3a, despite the fact that the butterfly hysteresis passes through the origin 
of the coordinates in Fig. 2. 

In conclusion, we have shown Fe 2 to exhibit unusually slow magnetization re- 
laxation effects below 5 K. These can be ascribed to an energy barrier A = 2.40 
K resulting from the S = 5 dimer ground state of [Fe 2 F 9 ] 3 " and a negative ax- 
ial anisotropy. The relaxation partly occurs via a thermally-assisted quantum 
tunneling process. Compared to other SMMs, the energy exchange between 
the spin system and the phonons is extremely inefficient in Fe 2 . Therefore 
the SMM phenomenon is observed at temperatures comparable to the barrier 
height and far above the blocking temperature, which for Fe 2 is expected to 
lie in the mK regime. To our knowledge the Fe 2 dimer is the smallest spin 
cluster so far to exhibit slow relaxation. 
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Figure Captions 



Fig. 1 Temperature dependence of the magnetic susceptibility yT of a sin- 
gle crystal of (Et4N) 3 Fe2F 9 in the range 1.8-50 K (circles). The crystal was 
oriented with its c axis parallel (H || c) and perpendicular (H _L c) to the 
magnetic field of 1 kG. The solid lines are calculated using eq. 1 with the 
parameters J = —2.23 K, D — —0.215 K, and g = 2.00. The inset shows the 
[Fe 2 F 9 ] 3- dimer molecule and the magnetic field direction for both H \\ c and 
H _L c orientations. 

Fig. 2 Positive wings of the butterfly hystereses of the magnetization in H \\ c 
orientation at 1.8 K for sweeping rates of T ^3 G/s (*) and T ml G/s (o). 
The 1.8 K magnetization in H _L c is shown in the inset. The solid lines are 
the adiabatic curves, calculated with J = —2.23 K, D — —0.215 K, and g = 
2.00 derived from susceptibility measurements. 

Fig. 3 a) Relaxation of the magnetization at fixed H z and T as indicated, 
after saturation at 30 kG and quick (T ps 150 G/s) sweep of the field to the 
indicated H z value. The solid lines are fits to a single exponential law. Note the 
different scale for the H z = measurements, b) Dependence of the relaxation 
time r from H z — with errors (3a). The dotted line is a guide to the eyes. 

Fig. 4 Magnetization data for H \\ c orientation at 1.8 K obtained as follows: 
After saturation at +30 kG, the field was quickly (T m 150 G/s) reduced to 
the given H z values including negative ones, and the magnetization was mea- 
sured within 50 s after reaching H z . The solid line is the calculated adiabatic 
magnetization curve at 1.8 K, see also Fig. 2. 
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Fig. 5 Schematic plot of the potential energy vs magnetization direction for 
a S = 5 ground state split by a negative axial zero-field splitting at H z = 
0. The potential energy barrier is A = (\D S=5 \)S 2 = 2.40 K. The full and 
broken arrows denote the thermal relaxation and possible tunneling pathways, 
respectively 
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